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We predict that the Josephson current in a clean s-wave 
superconductor-normal metal-d-wave superconductor junc- 
tion is periodic in superconducting phase difference ip with pe- 
riod TT instead of 2ti. The frequency of non-stationary Joseph- 
son effect is correspondingly 2ujj = AeV. The effect is due to 
coexistence in the normal layer of current carrying Andreev 
levels with phase differences ip and ip + n. 



The Josephson effect in superconductor- normal metal 
-superconductor (SNS) junctions differs from the effect 
in tunneling junctions in several important aspects [Q. 
In particular, in a long clean SNS junction at T = the 
Josephson current, Ij{ip), is a sawtooth function of super- 
conducting phase difference if g-Q] (solid line in Figjlp,) , 
in contrast to Ic sin ip dependence in tunneling junctions 
. The reason of this behaviour is that the Josephson 
current is transferred through the normal layer by current 
carrying Andreev levels, formed due to subgap Andreev 
reflections of electrons and holes from the non-diagonal 
pairing potential of superconductors |^,||. Positions of 
these levels, and therefore the current, depend on (p. 

In this Letter we consider the case when one of the 
superconductors has d-wave pairing symmetry: s-wave 
superconductor - normal metal -d-wave superconductor 
("SND") junction (FigHa). There is a strong evidence 
that such pairing is realized in superconducting cuprates, 
based particularly on phase sensitive measurements P,n| . 
The experimental realization of SND junction is feasible, 
and could provide additional information about super- 
conductivity in cuprates. Both tunneling and SNS junc- 
tions between conventional and d-wave superconductors 
are now being actively investigated |^. 

We predict that SND junctions will combine proper- 
ties of conventional and so called 7r-junction |9[[lCt] with 
negative critical current [|ll|(Fig.^, dashed line). As a 
result, they exhibit half-periodic Josephson effect, with 
/((/)) being 7r-periodic function of phase (instead of 2tt), 
which can be detected, e.g., by doubling of frequency of 
non-stationary Josephson effect. 

The physical picture of the effect is as follows. In d- 
wave superconductors the order parameter can have ei- 
ther sign depending on the momentum direction on the 
Fermi surface, which can be treated as an additional, in- 
trinsic phase TT |12]. 
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FIG. 1. (a) Current-phase characteristics of a clean 
SNS junction at zero temperature: conventional ("zero"-) 
junction (solid line); vr -junction (dashed line). (b) Cur- 
rent-phase characteristics of half-periodic Josephson current: 
Z = 0,T = (solid line); Z = 2^/^ ~ l,T = (dashed line); 
Z = 0,L/It = 0.2 (dash-dotted line); Z = 0, L/Zt = 0.5 
(dotted line). Disbalance function Z is defined in Eq.(9). (c) 
Current-phase characteristics of half-periodic Josephson cur- 
rent at finite normal scattering: Z — 0; L/U = 0.05 (dashed 
line); 0.1 (solid line). 

Therefore the Andreev levels in SND junction depend 
(instead of ip>) on effective phase difference 



(^(•^ii) ^ 



<p + TT d{~^q), 



(1) 



where §{x) is the Heaviside step function. Unit vector 
q gives the direction of the wave vector of transmitted 
state in d-wave superconductor, q, with components k|| 
(parallel to the interface) and qz{E) (FigJ3b). There are 
thus two sets of Andreev levels in the normal layer: with 
(p — ip (Ag positive) and (p = ip + t: (A^ negative). 
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FIG. 2. (a) s-Wave superconductor - normal metal - d-wave 
superconductor junction ("SND"-junction). We can choose 
real order parameter in d-wave superconductor, (b) Zero- 
and IT- Andreev levels in a clean SND junction. 

The Josephson current through 7r-levels will be a 2ti- 
periodic function of phase, shifted by tt with respect to 
the current carried by zero levels. The resulting current 
is TT-periodic in phase (Fig.|l|b). 

From the Josephson relation dtf = 2eV we see that 
the Josephson frequency doubles: a),/ — 2ujj — AeV. 

In our further analysis we assume that the following 
conditions are satisfied: 



max(^o) ■^L<s:lTje,k 



(2) 



where L is the width of the normal metal layer, max(^o) 
the larger of the correlation length in either supercon- 
ductor, li,l^ are elastic and inelastic scattering length 
respectively, and It — vp / {2'KkBT) is the normal metal 
coherence length in the clean limit. The latter is the 
characteristic coupling length in SNS junction ^^^ , and 
can be much larger than ^g (especially in case of cuprates 
[p^). In silver, e.g., at IK It — 1.67-10"'' cm, compared 
to ^0 ~ 10^^ cm typical for conventional superconduc- 
tors. 

Andreev levels in the normal layer are obtained by 
solving Bogoliubov-de Gennes equations for the two- 
component wave function in both superconductors and 
normal layer, separately for each k||-mode, and match- 
ing the wave functions at the interfaces |1J] ; for the SND 
junction we have: 



n 



A*W 



Here H 



A(z) 

-(n- 




= E 



u{z) 
v{z) 



2^V^ is the one-particle Hamiltonian; A(z) 
is the non-diagonal potential (which is zero in normal 



layer, A'^^^e*'^ in s-wave superconductor, and AV 



(2) 



in d- 

wave superconductor; we choose real A*^^\ A^ '). This 
standard approximation is justified by the condition ^q <C 
L [Q. We assume for simplicity the same value of kp 
in all three regions. The corrections due to differences in 
kp and finite normal scattering on the interface will be 
considered elsewhere. 

In the normal layer the normal component of the mo- 
mentum of an electron (hole) with energy E and tangen- 

HE) = ^ 



tial momentum kn is kl 



— \ I Kp 



/fcjj ± 2'mE. 



In 



the superconductor, for a subgap quasiparticle, E < |A|, 
it transforms into 



jt{E) = kM^ 



l±i 



lAI 



1- 



\ A;,(0)2/2mY |A| 



E^ 

AP- 



(4) 



If |A| ^ ^^i then ^qf « fcz(O), i.e. the quasiparticle 
momentum does not change direction in the supercon- 
ductor: 5 « fc, k = (kii , kz)- Therefore the momentum of 
the quasiparticle in the normal layer determines the value 
A j, of the d-wave order parameter entering Bogoliubov-de 
Gennes equations (y), and the effective phase difference 
for given Andreev level (Fig.||b). 

Under the above conditions the energies of low-lying 
Andreev levels [E <^ |A|) are given by 

{kl{E) ~ k'^{E)) L ± f>^^iO ^ TT{2n + 1); n = 0, ±1, ±2, . . 

(5) 

This is a direct generalization of the result by Kulik la] . 
The Josephson current through the contact can be cal- 
culated using the low-energy excitation spectrum deter- 
mined from (|) |,|l|, to yield at T = 0: 



Ijiv) = 



,-^ 2ev 



(k) 
Fz 



^^ nL 



n=\ 



(-ir 



-^ sva.nCp^'^i 



(6) 



Here ip'^'^^ is the effective phase shift, and v^l is the nor- 
mal to the interface component of Fermi velocity in a 

mode with tangential momentum kii — k,: I r 

•J\K\<~kp 



mv 



(k) 



= kp. The summation ^- 



dtZj^ (iK,y IS ex- 



tended over all allowed tangential modes fIB]; S is the 
area of the junction. 

It follows from (0) that the Josephson current in clean 
SND junction is indeed a sum of independent contribu- 
tions of zero- and 7r-modes: 



(3) 



Ij{ip) = /o (a'^+^F{ip) + a(-)F((^7r)) = (7) 

-_ I {F{^)+F{^7r) + Z {Fi^yFi^TT))) . (8) 





FIG. 3. On calculation of disbalance function at different 
orientations of d-wave superconductor with respect to the in- 
terface: white regions on the Fermi surface (in normal metal) 
correspond to zero modes, shadowed - to vr-modes, or vice 
versa, (a) Interface parallel to a? — y^-plane. (b) Interface 
normal to x^ — j/^-pIane; Q is the angle between the SN inter- 
face and the nodal direction in aj^-j/'^-plane. 



is the 27r-periodic 



Here F[x) = | ^- ^(-l)"+i5HH- 

sawtooth of unit amplitude; -^0 = 1; 'YIik, '^'fz ^^ the critical 
current of clean SNS junction at zero temperature [§,D; 

and weight factors a^^^ — ^- '^^fI/^k'^Fz^ where the 
zb-summation is extended over zero- and 7r-modes respec- 
tively. 

The ratio between the two contributions in (|7|) gener- 
ally depends on the orientation of the d-superconductor 
with respect to the interface, parametrized by some an- 
gles 9a, through the disbalance function. 



z{ea) 



v( + ) 



m 



-ry(-) 



(Oa), 



v(+K0a)+a(-)i9a)' 



\z{Oa)\ < 1; 



(9) 



Z = ±1 corresponds to purely zero- or n- junction pJlOl 
respectively; Z = when amplitudes of the contributions 
from zero- and 7r-levels are equal. 

Due to nonsinusoidal form of the current-phase char- 
acteristics, the current is not zero if Z{9a) = (as it 
would be in case of tunneling contact [ pO[ , with sin Lp in 
place of F{(p) ll^). On the contrary, the resulting cur- 
rent is TT-periodic in phase (Fig.np), and its critical value 
is halved. 

The simplest case is an SND junction where x^-j/^-plane 
of the d-wave superconductor is parallel to the interface 



(Fig.ga): then Z{9a) = by symmetry. For cuprates, 
this (OOl)-plane is also the easiest cleavage plane, which 
makes it the best candidate for experimental observation 
of the effect. 

Another possibility is to choose the interface normal 
to that plane. It is easy to see that in this case the 
disbalance function depends on a single angle 9 between 
the SN interface and the nodal direction in oi^-y^-plane 
(see Fig.|b), and Z(9) = ± {y/2cos[9 - f ] - l) ; (0 < 
9 < f). Since \Z{9)\ < (%/2- 1), the Josephson current 
in this case always contains a 7r-pcriodic component (see 

Fig.|i|b). 

It is worth noting that in case of "DND" junction 
between two d-wave superconductors the half-periodic 
Josephson effect is present as well, but it disappears at 
certain orientations of the crystals. If the junction inter- 
face is parallel to a;^— y^-plane of both superconductors, 
then Z{a) = ± (l - 4f ) . Here < a < f is the angle 
between nodal directions in the superconductors, and the 
effect is absent if a = or ^ . 

So far we did not consider the effects of nonmagnetic 
impurity scattering and finite temperature. 

At finite temperatures the condition E <Si jA^'^j will 
be met for almost all of the Fermi surface if 



knT. 



S<|A(i)|;max|A^')| 



(10) 



corrections to the current from the regions close to 
the nodal lines will be small by the same parameter, 
/csT/max |A. '\. This additional condition on tempera- 
ture is in fact less restrictive that the one following from 

i). 

Weak scattering can be taken into account via broad- 
ening of Andreev levels due to finite elastic scattering 
time, T = k/vp iQ. Using the approach of Q, we find 
for L/U.L/It < 1: 



Ij{v) 



2ev 



(k) CXD 



Fz 



ttL 



n=l 



(-1) 



n+1 



exp 



-2n-^ 



(«) 



{L /l!^^) sin mp'^^ 
smh{nL/l)p ) 



(11) 



,(S) 



,(«) 



It^^ are effec- 



(see FigJ|b,c). Here l'^' =h^,l'T 
five scattering and coherence length in mode k respec- 
tively. The effect survives finite temperature and weak 
elastic scattering in the normal layer. 

In the limit of strong elastic scattering, L ^ k, or 
at higher temperatures, L ^ It, the sawtooth current- 
phase dependence in long SNS junctions reduces to the 
sinusoidal one lia] , and the effect disappears due to can- 
cellation of contributions from zero- and 7r-modes. 

Note added: 
After this paper was submitted, the preprint [ pO[ ap- 
peared, where analogous conclusions were made in the 
case of short SND junction {L <^ niin(^o))- 
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